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structures  for  anny  applications.  In  these  events,  materials  are  subjected  to  large  strains, 
high  strain  rates  and  rapid  increase  in  temperature  due  to  thennoplastic  heating.  In  such 
severe  conditions,  overall  performance  is  determined  by  the  evolution  of  flow  stress, 
failure  initiation  and  propagation,  and  commonly  in  the  form  of  adiabatic  shear  banding. 

Some  of  2XXX  series  aluminum-copper  (Al-Cu)  alloys  are  recognized  for  their 
decent  ballistic  properties,  and  therefore  they  have  been  used  as  an  armor  material  for 
lightweight  U.S.  Army  vehicles.  Most  recently,  an  Al-Cu-Mg-Mn-Ag  alloy  labeled  as  A1 
2139-T8  has  been  developed  and  is  evaluated  by  the  U.S.  Anny  Research  Labs  because 
of  its  better  ballistic  properties  and  higher  strength  than  its  predecessors.  The  underlying 
microstructure  is  believed  to  be  the  key  element  for  this  superior  performance. 

The  goal  of  this  study  is  to  explore  the  effect  of  composition  and  microstructural 
features  on  overall  dynamic  material  behavior  by  examining  mechanical  and  deformation 
behavior  of  different  Al-Cu  material  systems.  Starting  from  the  pure  single  crystal  and 
polycrystalline  A1  structures,  and  adding  a  different  element  to  chemical  composition  in 
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data. 
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1 


CHAPTER  1 
INTRODUCTION 


1.1  Motivation  and  Scope 

Aluminum-copper  (Al-Cu)  alloys  are  broadly  used  in  military  and  aerospace 
applications  due  to  their  relatively  high  strength,  respectable  ballistic  performance  and 
lower  density  value  compared  to  other  armor  materials.  Military  vehicles  shown  in  Figure 
1 . 1  can  be  exposed  to  high  strain  rate  and  potentially  localized  deformation  accompanied 
with  rapid  increase  in  temperature,  and  therefore  armors  made  of  Al-Cu  alloys  are 
expected  to  endure  these  extreme  conditions. 


Figure  1.1.  Ml 09  howitzer  (left)  and  Marine  Corps.  Expeditionary  Fighting  Vehicle 
(right)  are  two  examples  of  military  vehicles  utilizing  aluminum-copper  based  armors 
which  are  subjected  to  harsh  conditions. 

2139  Al.  alloy  is  recognized  for  its  superior  ballistic  performance  compared  to 
other  Al-Cu  based  alloys  (Cho  and  Bes  2006)  (Placzankis  and  Charleton  2009).  In  their 
paper,  Placzankis  and  Charleton  qualitatively  showed  this  distinguished  ballistic 
performance  of  2139  alloy  compared  to  other  widely  used  Al.  alloys  which  have  decent 
ballistic  properties  (Fig. 2).  Although  this  notable  performance  in  damage  critical 
applications  is  well-known,  the  underlying  mechanism  for  this  behavior  is  not  well 
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understood  and  requires  a  methodological  study  as  a  function  of  both  composition  and 
microstructural  configuration.  This  study  is  an  attempt  to  reveal  the  effect  of  such 
functions  on  the  mechanical  response  of  2139  alloy. 


Minimum  Ultimate  Tensile  Strength  in  LT  Direction  (ksi) 


Figure  1.2.  The  experimental  test  results  of  different  aluminum  alloys  show  that  2139 
alloy  has  a  higher  minimum  ultimate  tensile  strength  and  a  better  ballistic 
performance.  (Placzankis  and  Charleton  2009) 

Current  thesis  study  is  composed  of  three  major  components  as  outlined  below:  (i) 
experimental  investigation  of  rate  and  temperature  dependent  material  behavior,  (ii) 
understanding  localized  deformation  behavior  in  the  form  of  ASB,  and  (iii)  constitutive 
modeling.  Considering  single  crystalline  and  polycrystalline  aluminum  as  base  materials, 
a  set  of  Al-Cu  alloy  systems  is  obtained  by  adding  a  different  element  in  each  stage 
(Table  1.1)  and  these  alloys  are  tested  in  quasi-static  and  dynamic  regime  (10~3s-1  to 
8000  s-1 )  at  different  temperatures  (ranging  from  room  temperature  to  320°)  to  discover 
their  relative  roles  on  strain  hardening,  strain  rate  sensitivity  as  well  as  their  coupling 
with  elevated  temperatures. 

In  high  strain  rate  deformation  events  such  as  ballistic  penetration  and  blast 
loading,  adiabatic  shear  banding  (ASB)  that  is  driven  by  the  onset  of  thennoplastic 
deformation  instabilities  emerges  as  a  the  key  deformation  and  failure  mechanism.  It  was 
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found  that  the  precipitates  may  have  an  important  effect  on  the  blocking  ASB  by  acting 
as  a  barrier  and  also  the  adiabatic  temperature  rise  is  less  than  the  matrix  containing  the 
precipitate  (Elkhodary,  et  al.  2009).  Therefore,  due  to  its  effect  on  ballistic  properties,  it 
is  crucial  to  investigate  the  ASB  phenomenon  in  different  alloys  to  understand  the 
influence  of  different  precipitates.  To  introduce  shear  banding  for  experimental  analysis, 
a  shear  compression  specimen  (SCS)  was  designed  (Rittel,  Lee  and  Ravichandran  2002) 
(Vural,  Molinari  and  Bhattacharyya  2010).  The  SCS  allows  the  material  to  experience 
large  shear  strains  along  a  narrow  gage  section  without  geometric  deformation 
instabilities  (such  as  necking  or  barreling)  and,  therefore,  is  considered  as  one  of  a  few 
specimen  geometries  suitable  for  exploring  ASB. 

Final  part  of  the  study  consists  of  the  constitutive  modeling  of  the  alloys  that  have 
been  experimentally  tested.  These  constitutive  models  act  as  a  link  between  experimental 
and  computational  work.  In  addition  to  modeling  the  experimental  data  with  purely 
phenomenological  Johnson-Cook  (JC)  model  (Johnson  and  Cook  1983),  dislocation 
mechanics  based  the  Mechanical  Threshold  Stress  (MTS)  (Follansbee  and  Kocks  1988) 
(Kocks,  Realistic  constitutive  relations  for  metal  plasticity  2001)  and  the  Zerilli- 
Amstrong  (ZA)  model  (Zerrilli  and  Armstrong  1987)  were  employed  to  construct  a 
dislocation  based  model  that  are  capable  of  capturing  experimental  data  with  the  least 
possible  deviation  over  a  wide  range  of  strain,  strain  rate  and  temperature. 

1.2  The  Alloys  and  Processing  Techniques 

This  section  describes  the  processing  techniques  for  the  alloys  that  are  tested  in 
this  study.  All  alloys  with  the  exception  of  single  crystalline  aluminum  are  prepared  by 
our  collaborative  team  in  material  science  division  led  by  Prof.  S.  Tin  (Table  1.1). 
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Mechanical  response  of  single  crystal  and  polycrystalline  aluminum  was 
investigated  by  a  former  team  member  (H.  Tran,  MS  Thesis,  IIT)  whose  results  are 
presented  in  this  thesis  study  to  provide  a  wider  perspective.  Single  crystal  aluminum  was 
tested  in  the  as-received  state  without  any  further  processing.  Polycrystalline  aluminum 
specimens  were  extracted  from  a  plate  which  was  processed  by  cold  rolling,  and  then 
annealing  to  refine  the  grain  structure.  All  alloys  were  melted  from  high  purity 
components  by  induction  heating,  and  casted  into  a  steel  mold.  Alloys  containing  Mg 
element  were  melted  in  an  argon  atmosphere  to  minimize  oxidation.  Precipitation 
hardenable  alloys  additionally  aged  at  a  low  temperature  of  320F  to  peak  hardness.  Each 
alloy  were  hot  rolled  and  annealed  to  produce  controlled  grain  refinement,  and 
recrystallization  kinetics  of  the  various  alloys  was  tracked  (Snyder  2012). 


Table  1.1.  A  summary  of  the  tested  alloys  with  their  individual  strengthening 
mechanisms  and  their  processing  procedure 


Alloy  (values  in  wt.%) 

Strengthening  Mechanism 

Processed  State 

[100]  and  [111]  Single 
Crystals 

Inherent  crystallography 

Bridgman  method, 

6N  purity 

Polycrystalline  A1 

Grain  boundaries 

Rolled  and 
recrystallized  plate 

Al-O.lCu 

Solid  solution 

Rolled  and 
recrystallized  plate 

Al-4.5Cu 

Heterogeneous  precipitation  of 
0’ 

Peak  Aged  at  160°C 
(320°F) 

Al-4.5Cu-0.5Mg 

More  uniform,  fine  dispersion 
of  0’ 

Peak  Aged  at  160°C 
(320°F) 

Al-4.5Cu-0.5Mg-0.3Mn 

0’,  A1  Cu  Mn  dispersoids 

Peak  Aged  at  160°C 
(320°F) 

Al-4.5Cu-0.5Mg-0.3Mn- 

0-3Ag 

All  above,  plus  Q  precipitation 

Peak  Aged  at  160°C 
(320°F) 
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1.3  Experimental  Procedure 

1.3.1  Quasi-Static  Experiments.  To  introduce  strain  rates  ranging  from  10-3  tolO0  s'1 
,  a  servo  hydraulic  test  frame  (with  load  capacity  of  535kN)(Fig.  1.3)  with  MTS  Flex  Test 
SE  Controller  with  high  temperature  MTS  extensometers.  For  testing  materials  in  this 
range,  cylindrical  specimens  machined  by  MMAE  Machine  Shop  with  an  L/D  ratio  of  1 .2 
were  used.  Due  to  small  size  of  the  specimens,  the  tested  specimens  were  sandwiched 
between  two  alumina  bars  whose  diameter  is  much  greater  than  the  specimen’s  diameter. 
Similarly  because  of  the  specimen  size  considerations,  extensometer  ends  were  attached 
to  alumina  bars  and  the  displacement  readings  were  taken  from  extensometer. 

High  temperature  quasi-static  tests  were  conducted  using  the  same  test  frame  and 
the  only  difference  in  the  setup  was  the  thennocouples  attached  to  specimen’s  and 
alumina  bars’  surface  via  high  temperature  cement.  A  ceramic  heater  attached  to  the  test 
frame  was  used  to  increase  temperature  of  the  specimen  to  the  desired  temperature.  The 
three  heating  units  inside  the  ceramic  heater  were  controlled  with  a  PID  controller  and  the 
three  thennocouples  supplied  the  temperature  readings  for  this  process.  An  average 
heating  sequence  takes  approximately  ten  minutes  to  reach  120  degrees  Celsius  and 
fifteen  minutes  to  reach  220  degrees  Celsius. 
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Figure  1.3.  The  MTS  quasi-static  testing  frame  is  shown  after  a  high  temperature  test. 

Furnace  door  is  open  and  the  thermocouples  are  taken  out 

1.3.2  Dynamic  Experiments.  The  Split  Hopkinson  Pressure  Bar  (SHPB)  setup  was 
used  to  conduct  high  strain  rate  compression  tests  on  Al-Cu  materials  and  to  obtain  the 
stress-strain  relationship  based  on  stress  wave  data  recorded  from  strain  gages  located  on 
input  and  output  bars.  Figure  1.4  shows  the  schematic  of  the  SFIPB  setup. 

A  striker  bar  launched  from  an  air  pressurized  gun  impacts  the  input  bar  and 
generates  a  compressive  stress  wave  (incident  wave)  that  propagates  along  the  input  bar. 
When  the  incident  wave  hits  the  interface  between  the  specimen  and  the  input  bar,  a 
portion  of  this  wave  reflects  as  a  tensile  wave  (reflected  wave),  and  the  other  compressive 
component  is  transmitted  through  specimen  to  output  bar,  thereby  deforming  the 
specimen  during  this  high-speed  process. 
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Strain  Gages 


Figure  1.4.  Schematic  of  the  SHPB  setup 


Each  of  these  stress  waves  (incident,  transmitted,  and  reflected  waves  (Fig.  1.5)) 
are  recorded  through  a  high  speed  digital  oscilloscope  by  using  the  signals  from  strain 
gages  that  are  attached  to  the  input  and  output  bars.  The  stress  waves  are  essential  in 
determining  dynamic  stress-strain  relationship  of  the  specimen  being  tested. 


oj- 


Figure  1.5.  Schematic  of  the  incident,  transmitted,  and  reflected  waves 
These  waves  are  essentially  longitudinal  elastic  waves  which  correspond  to  the 
motion  of  the  particles  back  and  forth  along  the  direction  of  wave  propagation  such  that 
the  particle  velocity  is  parallel  to  the  wave  velocity.  If  the  wave  is  compressive,  they  both 
have  the  same  sense;  if  it  is  tensile,  they  have  opposite  sense  (Meyers  1994). 
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Figure  1.6.  Propagation  of  an  elastic  compressive  wave  in  a  thin  cylindrical  bar 
Suppose  an  elastic  compressive  wave  is  travelling  through  a  bar’s  cross  section 
with  the  velocity  of  c0  in  the  longitudinal  direction,  as  shown  in  Fig.  1.6.  According  to 
the  Newton’s  second  law: 


d2u  (  doYY 

A.pSx.w  =  A[axx  +  Sx.— 


(7y 


(1.1) 


where  A  is  the  cross-sectional  area  of  the  bar,  p  is  the  density  of  the  bar,  Sx  is  an 
infinitesimal  length  of  the  bar  through  which  the  wave  is  passing,  u  is  the  particle 
displacement,  and  oxx  is  the  stress  at  the  left  side  of  the  cross-section.  After 
simplification: 


d 2  u  d  <rxx 

P'  dt2  dx 

The  one  dimensional  (longitudinal  direction)  Hooke’s  Law  is: 


(1.2) 


®xx  ~  E0.  £ 


du 

oxx^E^.— 


(1.3) 


where  E0  is  the  elastic  modulus  of  the  bar  and  the  £  is  the  elastic  strain  in  longitudinal 


direction.  Substituting  Eq.  1.3  into  the  Eq.  1.2  gives: 


The  wave  speed  in  the  bar  is  related  to  the  bar’s  density  and  elastic  modulus  as: 


(1.5) 

After  substituting  Eq.  1.5  into  Eq.  1.4,  one  dimensional  elastic  wave  equation  in  a  thin 
bar  can  be  obtained: 


d2u  2  d2u 
dt2  C°"dx2 


(1.6) 


Before  moving  on  to  the  SHPB  equations,  it  would  be  worthwhile  to  study  the  behavior 
of  the  wave  at  the  interfaces.  D’Alembert’s  solution  to  the  elastic  wave  equation  is: 

u  —  f(c0t  —  x)  +  F(c0t  +  x)  (1.7) 

where  u  is  the  displacement,  f  is  the  wave  travelling  in  the  direction  of  increasing  x,  F  is 
the  wave  travelling  in  the  direction  of  decreasing  x.  Considering  a  wave  travelling  in  the 
direction  of  decreasing  x: 


u  —  F(c0t  +  x) 

(1.8) 

du 

dF 

d(c0t  +  x ) 

—  r 

(1.9) 

dx 

d(c0t  +  x) 

dx 

du 

dF 

d(c0t  +  x) 

(1.10) 

dt 

d(c0t  +  x) 

St  =Fc“ 

The  last  two  equations  show  that  partial  differentiation  of  the  displacement  with  respect 
to  the  distance  and  time  can  be  related  to  each  other  as: 
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du  1  du 
dx  c0  dt 


(1.11) 


Recalling  that  strain  e  —  particle  velocity  v  —  and  substituting  them  into  Eq.  1.11: 


£  —  —  v 
co 

a  1 


(1.12) 


Finally,  when  Eq.  1.12  is  substituted  into  Eq.  1.5: 

a  —  pc0v  (1.11) 

When  an  incident  wave  of  magnitude  ol  propagating  along  bar  A,  resulting  in  a 
particle  velocity  Vj  reaches  the  interface  between  bars  A  and  B,  part  of  the  wave  with 
magnitude  oT  is  transmitted  to  bar  B,  and  the  other  part  with  magnitude  oR  is  reflected 
back  to  bar  A. 


A  V' 

1/ 

VT  B 

Vr  < 

Figure  1.7.  A  longitudinal  wave  reflects  at  the  interface  of  two  bars 


The  force  equilibrium  between  bars  is: 

A  a  (°7  +  ok)  —  Abot  (1.12) 

where  Aa  and  AB  are  the  cross-sectional  areas  of  the  two  bars.  Due  to  velocity  continuity: 

vi  +vR  —  vT  (1.13) 

Particle  velocity  can  be  linked  to  the  stress  in  the  bar  according  to  the  Eq.  1.11: 
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0/ 

V!  =  - 

Paca 


°R 

Paca 


(1.14) 


(7  t 

VT  =  - 

Pbcb 

where  p  and  c  are  the  density  and  the  wave  propagation  speed,  respectively.  The  negative 
sign  for  vR  is  that  a  positive  stress  causes  a  negative  particle  velocity  upon  reflection  and 
vice  versa  (Meyers  1994).  When  Eq.  1.14  is  substituted  to  the  Eq.  1.13: 

0/  <7r  O'f 

- = -  (1.15) 

Paca  Paca  Pbcb 

From  the  equations  1.12  and  1.13: 

or  _  2ABpBcB 
°i  ABpBcB+  AaPaca 

(1.16) 

aR  _  AbPbcb~AaPaca 
ai  Abpbcb+Aapaca 

where  Ape  is  called  the  impedance  of  the  bar. 

When  the  impedance  of  the  bar  A  is  greater  than  that  of  bar  B,  the  reflected  wave 

is  of  the  opposite  sign  to  the  incident  wave  and  when  it  is  less  than  the  impedance  of  bar 

B,  the  reflected  wave  is  of  the  same  sign  as  the  incident  wave.  However,  no  matter  the 

impedances  of  the  bars,  the  transmitted  wave  is  always  of  the  same  sign  as  the  incident 

wave. 


Considering  bar  B  as  a  free  end  with  zero  impedance  and  elastic  modulus,  the 


equations  become: 


(J  rr 

—  =  0 
0/ 
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0r 

°i 


-1 


(1.17) 


So  if  a  compressive  wave  is  propagating  to  a  free  end,  there  will  be  no  transmitted  wave 
and  the  reflected  wave  is  in  tension  with  the  same  magnitude  as  the  compressive  wave. 
Another  special  case  arises  when  bar  B  is  considered  as  a  rigid  body  with  infinite  elastic 
modulus: 


C 1 j1 


2 


0/ 


1 


(1.18) 


Thus  if  a  compressive  wave  is  propagating  to  a  rigid  body,  the  reflected  wave  is  also  a 
compression  wave  with  the  same  magnitude  as  the  incident  wave.  Last  special  case  is 
when  the  two  bars  have  the  same  impedance: 

( 'J'T 

—  =  1 
°I 

(1.19) 

Or 

—  =  0 
O'/ 


As  seen  from  these  equations,  when  the  two  bars  have  the  same  impedance,  there  will  be 
no  reflected  wave  and  the  transmitted  wave  will  have  the  same  sign  and  magnitude  with 
the  incident  wave. 

SHPB  setup  consists  of  a  compressive  elastic  stress  wave  propagating  towards  a 
sandwiched  specimen  whose  impedance  is  smaller  than  that  of  the  bars  (Fig.  1.8). 
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Strain  gage 


Strain  gage 


□ 

jp®mjs 

Figure  1.8.  Wave  propagation  in  SHPB  setup 


As  the  compressive  stress  wave,  the  incident  wave,  reaches  the  incident  bar- 
specimen  interface,  part  of  it  is  reflected  as  tensile  wave,  called  the  reflected  wave, 
because  the  bars  have  higher  impedance  than  the  specimen.  The  other  part  of  the  incident 
wave  will  be  transmitted  into  the  specimen  as  a  compressive  wave  and  when  this  wave 
reaches  the  specimen-transmitted  bar  interface,  part  of  it  will  be  reflected  back  into  the 
specimen  as  a  compressive  wave  due  to  specimen’s  smaller  impedance.  The  rest  of  the 
wave  will  enter  the  transmitted  bar  as  a  compressive  wave.  Strain  gage  attached  to  the 
incident  bar  will  record  the  incident  wave  and  the  reflected  wave,  whereas  the  other  strain 
gage  records  the  transmitted  wave.  D’Alembert’s  solution  to  the  wave  propagation  in 
SHPB  test  (Fig  1.9)  is: 


X 

>u2 

°R 

aT 

- ^  ^ 

>  77- 

^v2 


Figure  1.9.  One  dimensional  wave  analysis  in  SHPB  test 


U1  ~  f(X  ~  C0t)  +  FX  +  C0t  —  li;  +  uR 

u2  =  g(x  -  c0t)  =  uT 


(1.20) 


When  the  last  equation  is  differentiated  with  respect  to  x: 
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£l=  f  +F'  =  Sj+  £R 
e2=g'  =  eT 

Differentiating  Eq.  1.20  with  respect  to  time: 

Vi  =  Cotf'+F')  =  C0(£;  +£r) 

V2  =  -C0  g'  =  Cq£t 


(1.21) 


(1.22) 


As  the  specimen  is  fairly  shorter  compared  to  the  bars,  it  is  practical  to  neglect  the  stress 
wave  propagation  in  the  specimen  (Gray  2000).  Then  the  engineering  strain  rate  of  the 
specimen  is: 

.  -  v2 


c  =  y~  (-£/  +  cR  +  er) 

^ S 


(1.23) 


£  = 


hcE, 


'  (“0/  +  +  0r) 


5^0 


where  Ls  is  the  original  specimen  length,  vx  and  v2  are  the  particle  velocities.  The 
engineering  strain  can  be  found  by  integrating  the  strain  rate: 


-r* 

co  r 1 
LSE0  J 


(1.24) 


(  —  (7;  +  Or  +  <JT)dt 


The  engineering  stress  can  be  calculated  by  assuming  that  the  force  acting  on  specimen  is 
the  average  of  the  two  forces  acting  on  the  two  interfaces: 
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All  the  derived  equations  are  based  on  the  one  dimensional  elastic  wave  propagation  in 
the  bars,  for  this  assumption  to  be  valid;  several  conditions  should  be  satisfied  (Gama, 
Lopatnikov  and  Jr.  2004): 

a)  The  bars  are  homogenous  and  isotropic. 

b)  The  bars  are  under  elastic  defonnation. 


c)  The  bars  are  free  of  dispersion. 
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1.4  Adiabatic  Shear  Banding 

Adiabatic  shear  banding  (ASB)  is  the  major  deformation  mode  in  high-strain-rate 
deformation  events,  such  as  ballistic  penetration  and  blast  loading,  it  is  particularly 
important  in  dynamic  response  and  failure  behavior  parameters  on  the  onset  of 
thennoplastic  deformation  instabilities  (Fig.  1.10)  (Meyers  1994). 


Figure  1.10.  Defeat  of  armor  by  plugging  is  an  example  to  formation  of  shear  bands  in 
dynamic  deformation  events. 


Adiabatic  shear  band  formation  involves  a  narrow  localized  defonnation  band 
where  very  high  shear  strain  takes  place  during  dynamic  deformation  processes  such  that 
there  is  not  enough  time  for  thermoplastic  heating  to  dissipate  away  from  localized 
deformation  site.  This  usually  follows  the  situations  where  a  bullet  or  a  projectile  has 
impacted  a  target  material  or  during  the  manufacturing  process  where  chipping,  forging 
and  upsetting  can  create  unwanted  adiabatic  shear  bands  (Wright  2002).  For  aluminum 
materials  adiabatic  shear  bands  even  occurred  in  cylindrical  compression  specimens 
(Wulf  1978). 

A  specialized  specimen  geometry,  called  shear-compression  specimen  (SCS),  has 
been  developed  to  study  formation  of  ASB  (Rittel,  Lee  and  Ravichandran  2002)  (Rittel, 
Ravichandran  and  Lee  2002).  A  cylindrical  specimen  geometry  (Fig.  1.11),  with  two 
diametrically  symmetric  slots  oriented  with  a  specific  degree  to  the  axis  of  specimen  was 
developed  to  stimulate  uniform  shear  strain  and  stress  distribution  along  the  gage  section. 
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More  recently,  the  specimen  geometry  was  improved  by  discovering  that  a  slot  angle  of 
35.26  degrees  promotes  the  desired  conditions  in  deformation  of  the  gage  section  (Vural, 
Molinari  and  Bhattacharyya  2010). 


Figure  1.11.  Schematic  of  SCS  geometry 


Equivalent  stress  and  strain  can  be  computed  from  the  measured  load  and 
displacement  as: 


_  1  8 

£eq  ~  VJi 


(1.27) 


oeq  =  kx  exp(-k2seq)  — 

where  8  and  are  measured  displacement  and  force,  respectively.  The  constants  /q  and 
k2  are  found  by  comparing  cylindrical  non-slotted  specimen  test  with  shear-compression 


specimen  test,  both  of  which  are  carried  out  at  similar  quasi-static  strain  rates. 
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1.5  Modified  Johnson-Cook  Model 

The  Johnson-Cook  (JC)  phenomenological  plastic  flow  stress  model  is  one  of  the 
most  widely  known  temperature  and  strain  rate  dependent  model  and  it  is  commonly  used 
in  computational  codes  (Johnson  and  Cook  1983).  It  has  the  following  relation  for  the 
flow  stress: 


a  =  K  +Bi£P)n] 


1  +  Cln 


[1  -  (T*y] 


(1.28) 


where  ep  is  the  equivalent  plastic  strain,  £  is  the  strain  rate,  £q  is  a  user  defined  reference 
strain  rate,  and  er0  ,  B,  n,  C,  p  are  material  constants.  The  nonnalized  temperature,  T* ,  is 
defined  as: 


T-Tn 


(1.29) 


Tm  To 

where  T0  is  a  reference  temperature  and  Tm  is  a  reference  melting  temperature  in  K.  The 
success  of  JC  model  comes  from  multiplicative  decomposition  of  strain  hardening,  strain 
rate  hardening  and  thermal  softening  effects  in  three  separate  terms.  Even  if  JC  model  is 
very  flexible  in  modeling  a  variety  of  materials,  it  fails  to  capture  the  following  behaviors 
of  aluminum-cupper  alloys; 

•  The  amount  of  strain  hardening  (B)  decreases  faster  with  increasing  temperature 
than  the  predicted  rate  in  JC  model. 

•  The  strain  rate  sensitivity  in  quasi-static  regime  increases  with  temperature, 


whereas  JC  model  employs  temperature  independent  strain  rate  sensitivity  in 
quasi-static  regime. 
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•  Regardless  of  the  temperature,  there  is  a  noteworthy  dissimilarity  in  strain  rate 
sensitivities  of  quasi-static  and  dynamic  regime.  However  JC  model  predicts 
identical  strain  rate  sensitivity  in  both  regimes. 

A  modified  Johnson-Cook  Model  (MJC)  was  developed  to  address  these  concerns  in 
JC  model  (Vural  and  Caro  2009).  The  strain  hardening  tenn  (B)  was  coupled  with 
temperature  as: 

B  —  B0[l  —  (T*)p]  (1.30) 

To  resolve  remaining  issues  in  JC  model,  the  strain  rate  sensitivity  term  (C) 
redefined  as: 


C  =  CiCr;)4*  +  C2  H(£,£t,k)  (1.31) 

where 


T~T0 

Tr~T0 


(1.32) 

11  /  £\ 

H(s,£t,  k)  —  —  +  -tanh  kin  — 

2  2  V  / 

Here,  £t  represents  the  transition  strain  rate  which  separates  the  quasi-static  and  dynamic 
regime.  H(£,£t,k)  function  is  a  Heaviside  step  function  that  gives  1  when  strain  rate  is 
greater  than  the  transition  rate  and  gives  0  when  strain  rate  is  smaller.  Moreover  by 
supplying  continuous  variation  around  the  transition  strain  rate,  this  step  function  ensures 
a  smooth  transition  between  the  two  different  strain  rate  regimes.  Thanks  to  this  step 
function,  only  C1  term  dominates  the  quasi-static  regime  and  also  since  it  is  coupled  with 
temperature,  it  makes  certain  that  there  is  increased  rate  sensitivity  with  increasing 
temperature  in  quasi-static  regime.  Furthermore,  C2  can  be  chosen  to  be  much  greater 
than  the  Cx  and  thus  dynamic  regime  can  be  controlled  merely  by  C2  tenn.  Lastly,  due  to 
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continuity  in  the  flow  stress  at  the  end  of  quasi-static  regime  and  at  the  beginning  of 
dynamic  regime  requires  that 
Gqs  d  ^ 


£n  —  £n 


C2 

vC1(Tr)q+C2 


-1  H(ep,i,k) 


(1.33) 


where  £01  is  the  reference  strain  rate. 


1.6  Mechanical  Threshold  Model 


Mechanical  Threshold  Stress  (MTS)  model  is  a  dislocation  based  material  model 
which  incorporates  mechanical  threshold  stress  as  a  state  variable  (Follansbee  and  Kocks 
1988)  (Follansbee,  Kocks  and  Regazzoni  1985)  (Kocks,  Realistic  constitutive  relations 
for  metal  plasticity  2001).  .  To  comprehend  the  meaning  of  the  mechanical  threshold 
stress,  a  thought  experiment  can  be  carried  out  at  0  K.  If  an  increasing  amount  of  shear 
stress  is  applied  to  a  nearly  perfect  crystal  containing  the  necessary  defects  for 
deformation,  there  would  be  no  detectable  shear  strain  rate  until  a  =  f  (0).  At  this  shear 
stress,  the  shear  strain  rate  will  approach  an  inertially  limited  strain  rate,  and  this  material 
property  f(0)  is  called  as  the  mechanical  threshold  for  initiating  deformation  by  a 
specific  mechanism,  such  as  the  intrinsic  lattice  resistance  to  dislocation  glide  or 
resistance  to  dislocation  glide  caused  by  a  field  of  solute  atoms,  precipitate  particles. 
Therefore  this  idealization  can  be  summarized  at  0  K  as: 


y  =  0  for  a  <  f (0), 

Y  >  0  for  o  —  f(0). 

also  a  —  f  (0)  is  unattainable  under  quasi-static  conditions. 


(1.34) 


When  temperature  is  greater  than  zero  Kelvin,  this  behavior  will  change  in  two 
ways.  Since  the  resistance  to  deformation  is  governed  by  the  elastic  interaction  on  atomic 


21 


level,  temperature  dependence  of  the  elastic  constants. (Where  is  the  sentence  here?) 
Similarly,  with  increasing  temperature,  local  energy  barriers  can  be  overcome  at  lower 
shear  stress  with  the  help  of  the  local  temperature  fluctuations. 

The  Mechanical  Threshold  Stress  Model  uses  this  mechanical  threshold  stress 
(flow  stress  at  0  K)  as  a  state  parameter  to  calculate  the  flow  stress.  The  mechanical 
threshold  stress  is  separated  into  two  components, 

<7  =  da  +  at  (1.35) 

where  oa  represents  the  rate  independent  interactions  of  dislocations  with  long  range 
barriers  such  as  grain  boundaries  whereas  dt  represents  the  rate  dependent  interactions 
with  short  range  obstacles.  The  flow  stress  at  temperatures  greater  than  0  K,  the  flow 
stress  will  be: 


g  =  da  +  ot  =  da  +  s(e,  T)  dt  (1.36) 

If  the  rate  controlling  deformation  mechanism  is  thennal  activation,  then  the  contribution 
of  thermal  activation  energy  will  reduce  the  stress  required  for  a  dislocation  to  pass  a 
barrier.  Thus  the  factor  S  always  will  be  less  than  1 .  The  kinetics  of  thermally  activated 
glide  is  described  by  the  Arrhenius  expression  of  the  form: 

\AG(at,dt)] 


e  =  60  exp 


kT 


(1.37) 


AG  =  g0  gb: 


-w 


(1.38) 


For  the  free  energy  AG,  a  phenomenological  relation  is  selected: 

vP-|4 

>t} 

where  g0  is  the  normalized  activation  energy,  which  is  assumed  to  remain  constant  as 
long  as  obstacle  character  is  same,  g  is  the  shear  modulus,  b  is  the  magnitude  of  Burgers 


vector,  and  p  and  q  are  constants  that  characterize  the  obstacle  profile.  For  different 
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materials  and  for  different  heat  treatments  for  the  same  material  p  and  q  are  different.  If 
last  three  equations  are  rearranged,  the  equation  for  flow  stress  will  be: 


<*  =  °a  +  (#  “  #a)  \  1  ~ 


f 

1  - 

kTln  (V) 

1\ 

Q 

g0pb3 

V 

y 

(1.39) 


This  equation  gives  the  flow  stress  value  for  a  given  structure  (a),  thus  investigating  the 
evolution  of  the  structure  is  necessary  for  obtaining  stress  strain  curves. 

MTS  model  uses  Kocks-Mecking(KM)  model  (Kocks  1976)  for  predicting  the 
evolution  of  the  structure  with  the  strain.  This  approach  assumes  that  the  kinetics  of 
plastic  flow  is  determined  by  a  single  structure  parameter  S  representing  the  current 
structure.  At  a  given  structure,  flow  stress  can  be  represented  as: 

A  a  =  a(S,s,T)  (1.40) 

Evolution  of  the  structure  with  strain  £  at  a  given  temperature  and  strain  rate  is: 

=  f(S,s,T)  (1.41) 

as 

It  is  also  assumed  that  the  structure  parameter  S  evolves  towards  a  saturation  value  Ss. 
Thus,  the  flow  stress  o  evolves  towards  a  saturation  or  steady  state  of  value  os.  as  can  be 
expressed  as  a  function  of  temperature  and  strain  rate  as  os  —  os  (i,  T). 

MTS  model  chooses  the  mechanical  strength  of  obstacles  to  dislocation  glide,  a, 

as  the  structure  parameter  which  is  related  to  dislocation  density  as: 

i 

a  —  agbpz 


(1.42) 
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where  p  is  the  shear  modulus,  b  is  the  magnitude  of  the  Burgers  vector  and  a  is  a 
numerical  constant  of  order  unity.  The  KM  model  has  the  equation  for  the  evolution  of 
dislocation  density,  p,  as: 

dp  r 

—  =  /qp2  -  k2p  (1-43) 

de 

1 

The  ktp2  tenn  is  associated  with  the  athennal  ( k1  ^constant)  storage  of  moving 

dislocations  which  become  immobilized  after  having  travelled  a  distance  proportional  to 

1 

average  spacing  between  dislocations,  pT  k2p  is  the  dynamic  recovery  by  dislocation 
annihilation  and  cross  slip  of  the  dislocations.  Since  dynamic  recovery  is  thermally 
activated,  k2  is  a  function  of  temperature  and  strain  rate. 

By  combining  the  last  two  equations,  the  evolution  equation  for  the  structure  can 
be  found  as: 


(1.44) 


with 


9  = 


apbk-t 


CTc 


= af,b  (I) 


(1.45) 


According  to  the  KM  model  9  does  not  change  with  strain  rate  and  has  temperature 
dependence  through  p  term.  Also  according  to  this  equation,  d  evolves  towards  a  steady 
state  value,  ds  where  strain  hardening  is  zero.  The  MTS  model  did  not  employ  the  KM 
hardening  model  without  any  modifications.  First  to  have  a  better  fit  to  the  data,  it  has 


been  modified  to: 
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e  =  0O  [1  5-) 

L  V(7S  -  (7a/ 


(1.46) 


The  function  F  is  selected  to  fit  the  measured  data.  Another  deviation  from  the  original 
KM  hardening  model  is  that  the  athermal  hardening  coefficient  6{)  is  found  to  have  a 
strain  rate  dependency  from  the  experimental  data  in  FCC  materials.  Thus,  it  is  assumed 
have  a  form  such  as: 


90  =  A  +  Bln(s)  +  C  i  (1.47) 

The  KM  hardening  model  assumes  the  dislocations  move  a  distance,  d,  proportional  to 
the  average  distance  between  dislocations.  But  at  very  high  strain  rates,  the  dislocation 
immobilization  distance  is  simply  the  distance  that  a  dislocation  can  move  during  the 
imposed  time  duration  of  deformation.  Also,  this  increase  in  athermal  strain  rate  with 
strain  rate  may  be  the  reason  for  the  increased  strain  rate  sensitivity  in  dynamic  region  for 
fee  materials.  Finally,  an  Arrhenius  type  of  equation  is  proposed  for  saturation  stress: 


In 


(1.48) 


where  £s0,  <rs0,  A  are  constants  and  ds  is  the  saturation  threshold  stress  for  deformation  at 


OK. 


After  recent  modifications,  latest  fonn  of  the  MTS  model  is  (Banerjee  2007): 

g  =  Ga  +  (SjCXj  +  Seae )  ^ V ’  ^  (1.49) 

Mo 

where  Ga  is  the  athermal  component  of  stress,  gl  is  the  intrinsic  component  of  the  flow 
stress  ,  Ge  is  the  strain  hardening  component  of  the  flow  stress,  5)  and  Se  are  strain  rate 
and  temperature  dependent  scaling  factors,  and  /r0  is  the  shear  modulus  at  0  K  and 
ambient  pressure.  Parallel  to  the  original  MTS  model  (Eq.1.39),  scaling  factors  have 


Arrhenius  fonn: 
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i 


St 


Se 


f 

1  \ 

kTln 

Qi 

1  - 

V  £  J 

goikb3 

V 

y 

r 

1  ^ 

kTln  (^) 

Qe 

1  - 

goekb3 

■ 

V 


J 


(1.50) 


where  g oie  is  the  normalized  activation  energy,  ,  g  is  the  shear  modulus,  b  is  the 
magnitude  of  Burgers  vector,  £oi  e  are  the  reference  strain  rates,  and  pi  e  and  qLe  are 
constants  that  characterize  the  obstacle  profile. 

Hardening  is  kept  same  as  the  original  equation  with  the  addition  of  saturation 
hardening: 


^=eo[l-F(<re)]  +  01F(oe) 

90  —  a00  +  awlne  +  a20VZ  +  a30F 

91  —  a01  +  a^lns  +  a21V^  +  a31T 

tanh  (a  — ) 
ur  -  1  °sJ 

tanh  (a) 


(1.51) 


In 


where  0O  is  the  hardening  rate,  Oj  is  hardening  rate  at  saturation,  ai;-,  and  a  are  constants, 
£s0,  <Jes,  A  are  constants  at  zero  strain  hardening  rate  and  ae0s  is  the  saturation  threshold 
stress  for  deformation  at  0  K. 
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1.7  Zerilli-Armstrong  Model 

The  Zerilli-Armstrong  (ZA)  model  (  (Zerrilli  and  Armstrong  1987)  (Zerrilli  and 
Armstrong  1997)  (Zerilli  2004))  is  based  on  simplified  dislocation  mechanics. 
Dislocation  model  basis  for  plastic  shear  strain-rate  is: 

Y  —  m'bpv  (1.52) 


where  m’  is  a  tensor  orientation  factor,  p  is  the  dislocation  density,  b  is  the  Burgers 
vector  and  v  is  the  average  dislocation  velocity.  Average  dislocation  velocity  is  governed 
by  thermally  activated  processes  to  overcome  local  obstacles: 


v  —  v0  exp 


(1.53) 


where  v0  is  the  reference  dislocation  velocity,  G  is  the  Gibbs  free  energy  of  activation,  k 
is  the  Boltzmann  constant,  and  T  is  the  absolute  temperature.  G  can  be  expressed  as: 

rTth 

G  =  G0-  A*b  dx'th  (1.54) 

7o 

where  G0  is  the  reference  Gibbs  energy  at  T=0,  A*  is  the  area  of  activation,  and  r^is  the 
thennal  component  of  the  shear  stress.  Mean  value  of  A*b  can  be  expressed  as: 


<  Ab  >=  (  — 

^ th 


Tth 


A*b  drl 


th 


(1.54) 


Eq.  1.52  to  Eq.1.54  can  be  expressed  in  an  alternative  for,  keeping  in  mind  that  ath  = 


m  xth  and  e  —  — : 

Ln  m 


ath  =  Be  PT 


(1.55) 


where 
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B  = 


mG0 

A0b 


(1.56) 


The  B  term  is  the  threshold  stress  and  it  is  reduced  by  e  PT  factor  to  calculate  the  flow 
stress.  Using  the  expansion  ln(l  +  x)  ~  x  ,  the  /?  tenn  can  be  written  as: 

/?  =  /?0-/?1ln(f)  (1.57) 

Activation  area  (A)  is  a  measure  of  the  separation  distance,  d,  between  the  dislocations 
such  as: 


(1.58) 


Also  it  is  known  that  average  dislocation  separation  (d)  is  related  to  dislocation  density 
(P)  as: 


1 


p~ 


d2 


(1.59) 


Using  Eq.1.59,  activation  area  can  be  related  to  dislocation  density: 

1 


(1.60) 


Considering  the  relation  between  strain  and  dislocation  density  at  0  K: 

A~Te  (L61) 

Finally  activation  energy  can  be  described  with  following  relation: 

A  =  ^L  (1.62) 

y]e 

Thus,  for  the  fee  case,  thermal  component  of  the  stress  can  be  expressed  as: 
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ath  —  CiVeexpC— czT  +  c3Tlne )  (1.63) 

There  will  be  an  effect  of  the  grain  boundaries  on  the  stress  and  ZA  model  covers  this 
effect  by  adding  an  athermal  component  of  stress  to  the  model: 


a“  ^i+as 


(1.64) 


Where  k  is  the  microstructural  stress  intensity,  1  is  the  average  grain  diameter  and  os  is 
the  stress  increase  due  to  solutes.  Final  form  of  the  ZA  model  is: 


_ _  rv 

a  =  ctyje  exp(— c2T  +  c3Tlne )  +  —  +  <js 

V/ 


(1.65) 


Zerrilli  (Zerilli  2004)  modified  this  equation  to  include  both  Peierls  stress  type 
interactions  and  interaction  of  forest-dislocations  type  interactions: 

o  =  oa  + Be~PT  +  Bq£  2e~aT  (1.66) 

where 


a  =  a0  —  a1  ln(f) 

P=Po~Pi  ln<A) 


On 


(1.67) 
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CHAPTER  2 

EXPERIMENTAL  RESULTS 

2.1  Ultra-High-Purity  Single  Crystal  Aluminum 

As  a  first  step  in  this  study,  uniaxial  compression  tests  were  conducted  in  quasi¬ 
static  and  dynamic  regimes  at  room  temperature  in  ultra-high-purity  single  crystal 
aluminum  in  the  [111]  and  [001]  orientations  (Tran 201  l)(Fig.  3.1)(Fig  3.2). 


Figure  2.1.  On  left  force  going  into  [111]  orientation  and  on  right  going  into  [001] 
orientation  of  face-centered  cubic  structure.  [Provided  by  H.  Tran] 

As  figure  2.2  shows  that  the  [111]  orientation  single  crystal  has  twice  as  larger 
flow  stress  values  compared  to  [001]  orientation  single  crystal  at  a  given  strain  value. 
Also  both  orientations  show  increased  rate  sensitivity  in  dynamic  region.  They  show 
continuous  strain  hardening  without  saturation  at  large  strain  values.  Since  there  were 
limited  number  of  single  crystal  aluminum  specimens,  no  high  temperature  tests  were 
done,  thus  this  material  was  not  characterized  at  elevated  temperatures. 
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Single  Crystal  Al  111 


Single  Crystal  Al  001 


(111)  Single  Crystal  Al  Rate 
Sensitivity 


(001)  Single  Crystal  Al  Rate 
Sensitivity 
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Figure2.2.  The  material  response  of  the  two  single  crystals  at  quasi-static  and  dynamic 
regimes  can  be  seen  on  the  top  row  and  the  strain  rate  sensitivity  of  the  two  single 
crystals  is  presented  at  lower  row.  [Provided  by  H.  Tran] 


2.2  High-purity  Polycrystalline  Aluminum 

After  testing  the  single  crystalline  aluminum,  the  next  step  was  testing 
polycrystalline  aluminum  material.  As  polycrystalline  aluminum  contains  various  grains 
in  different  orientations,  the  material  strength  would  be  somewhere  in  between  the  two 
single  crystal  aluminums  (Fig.  2.3) 
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Figure  2.3.  Comparison  of  the  quasi-static  performance  of  single  crystalline  and 
polycrystalline  aluminum 


Tests  in  quasi-static  and  dynamic  regimes  were  conducted  at  room  temperature 
and  elevated  temperatures  (120°  C  and  220°  C)  (Fig.  3.5). 

2.3  Solid  Solution  Aluminum  Copper  (0.1%  Cu)  Alloy 

Next  alloy  is  obtained  by  adding  0.1%  (by  weight)  copper  to  the  polycrystalline 
aluminum,  which  resulted  in  an  aluminum  alloy  strengthened  by  solid  solution  of  copper 
atoms.  Flow  stress  values  are  very  similar  to  that  of  polycrystalline  aluminum  (Fig  3.4). 
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Figure  2.4.  In  quasi-static  regime  solid  solution  aluminum  copper  alloy  display  slightly 
higher  flow  stress  values  at  0. 1  strain  offset;  while  in  dynamic  regime  both  materials 
have  very  similar  flow  stress  values  (  Also  similar  trends  in  elevated  temperatures). 
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Figure  2.5.  Top  row  and  lower  left  graphs  show  the  material  response  of  polycrystalline 
aluminum  in  different  temperatures.  Lower  right  graph  demonstrates  the  strain  rate 
sensitivity  of  the  material  in  these  temperatures  at  0. 1  strain  value. 


To  fully  characterize  the  material,  tests  were  done  in  quasi-static  and  dynamic 
regimes  at  room  temperature  and  elevated  temperatures  (120°  C  and  220°  C)  (Fig.  3.6). 
Material  response  of  polycrystalline  aluminum  and  solid  solution  aluminum  copper  alloy 
are  almost  identical  and  they  show  increased  rate  sensitivity  at  the  dynamic  regime.  Both 
materials  also  have  increasing  strain  rate  sensitivity  with  temperature  in  the  quasi-static 
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regime.  Room  temperature  flow  stress  values  of  these  materials  is  much  closer  to  [001] 
single  crystalline  aluminum  than  the  [111]  orientation. 
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Figure  2.6.  Top  row  and  lower  left  graphs  show  the  material  response  of  solid  solution 
aluminum  copper  alloy  in  different  temperatures.  Lower  right  graph  demonstrates  the 
strain  rate  sensitivity  of  the  material  in  these  temperatures  at  0.1  strain  offset  value. 
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2.4  Al-4.5%Cu  Alloy 

With  the  addition  of  the  4.5%  Cu,  the  new  alloys  start  to  become  more 
microstructurally  complex  compared  to  the  previous  materials  which  do  not  contain  any 
precipitations.  Al-4.5%Cu  alloy  is  the  first  alloy  to  have  precipitations  which  are  in  the 
form  of  heterogeneous  precipitation  of  0’  phase.  Compared  to  previous  materials  which 
have  flow  stress  values  less  than  100  MPa,  flow  stresses  jump  to  400  MPa  with  the 
addition  of  precipitations.  At  this  time  only  room  temperature  data  is  available  for  this 
alloy. 
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Figure  2.7.  Room  temperature  experimental  data  for  Al-4.5%Cu  is  shown.  Since 
experiments  are  done  till  3000/s  strain  rate,  flow  stress  values  are  very  similar  in  all 
experiments. 
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2.5  Al-4.5%Cu-0.5%Mg  Alloy 

The  next  alloy  is  obtained  by  adding  0.5%Mg  (by  weight)  to  the  previous  alloy.  0' 
precipitations  that  were  previously  fonned  by  the  addition  of  the  4.5%Cu  becomes  more 
uniform  with  the  addition  of  Mg  and  the  new  alloy  comes  out  to  be  stronger  than  the  Al- 
4.5%Cu  alloy  (Fig.  2.8). 


AI-4.5%Cu  v.s.  AI-4.5%Cu-0.5%Mg 


Figure  2.8.  Comparison  of  Al-4.5%Cu  and  Al-4.5%Cu-0.5%Mg  at  the  same  strain  rate. 
Addition  of  the  Mg  generates  higher  flow  stress  values. 

At  this  point,  only  room  temperature  experimental  data  in  quasi-static  and 
dynamic  regimes,  and  220°  C  quasi-static  regime  experimental  data  is  available  (Fig. 
2.9). 

2.6  Al-4.5%Cu  0.5%Mg-0.3%Mn  Alloy 

Mn  element  (0.3%  by  weight)  was  added  in  this  stage  to  the  earlier  Al-4.5%Cu- 
0.5%Mg  alloy  which  resulted  in  Al2oCu2Mn3  dispersoids  and  0'  precipitations  to  exist  in 
this  latest  Al-4.5%Cu  0.5%Mg-0.3%Mn  alloy.  Flow  stress  is  slightly  higher  compared  to 
Al-4.5%Cu  0.5%Mg  alloy  (Fig  2.10). 
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Figure  2.9.  Current  progress  in  Al-4.5%Cu-0.5%Mg  alloy  is  presented.  Since  only  strain 
rates  up  to  3000/s  have  been  covered,  flow  stress  values  in  dynamic  regime  are  not 
very  different  from  the  quasi-static  flow  stress  values.  Only  quasi-static  data  is 
available  at  220°  C. 


Only  room  temperature  experimental  data  is  available  for  the  Al-4.5%Cu  0.5%Mg- 
0.3%Mn  alloy  in  the  quasi-static  and  the  dynamic  regime  (Fig.2.11). 
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AI-4.5%Cu  v.s.  AI-4.5%Cu-0.5%Mg 
v.s.  AI-4.5%Cu-0.5%Mg-0.3%Mn 


Figure  2.10.  Comparison  of  the  last  three  alloys  in  the  dynamic  regime.  Al-4.5%Cu- 
0.5%Mg-0.3%Mn  alloy  has  the  highest  flow  stress  values. 
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Figure  2.1 1.  Room  temperature  data  for  the  A-4.5%Cu-0.5%Mg-0.3%Mn  alloy  is 
presented. 
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2.7  Al-4.5%Cu  0.5%Mg-0.3%Mn  -0.3%Ag  Alloy 

The  most  complex  alloy  was  obtained  with  the  addition  of  Ag  element  which 
resulted  in  Q  precipitations,  AI20CU2M113  dispersoids  and  9'  precipitations.  Dynamic  and 
quasi-static  tests  were  conducted  at  room  temperature  and  220°C.  Also  0.001/s  quasi¬ 
static  tests  were  done  at  elevated  temperatures  to  capture  thermal  softening  behavior. 
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Figure  2.12.  Room  temperature  and  220°C  experimental  data  are  presented  at  the  upper 
row.  Lower  left  graph  show  the  thermal  softening  behavior.  Lower  right  graph  show 
all  the  experimental  data  in  0. 15  strain  offset. 
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Al-4.5%Cu  0.5%Mg-0.3%Mn  -0.3%Ag  alloy  has  the  highest  flow  stress  values 
among  the  alloys  that  were  tested  (Fig.  2.13). 


Figure  2.13.  Al-4.5%Cu  0.5%Mg-0.3%Mn  -0.3%Ag  alloy  has  the  highest  flow  stress 
values  among  the  other  complex  alloys.  However  after  0.3  strain,  Al-4.5%Cu 
0.5%Mg-0.3%Mn  reaches  the  same  flow  stress  value  with  the  Ag  alloy. 
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CHAPTER  3 

MODIFIED  ZERILLI-ARMSTRONG  MODEL 

3.1  Introduction 

Although  Modified  Johnson-Cook  model  excels  at  capturing  material  behavior, 

but  it  has  no  physical  roots.  MTS  model  seems  as  an  ideal  candidate  for  this  task,  but  it 

has  its  shortcomings.  Since  MTS  model  employs  threshold  stress  to  model  the  structure, 

first  threshold  stress  should  be  modeled.  To  model  evolution  of  threshold  stress  with 

strain  rate  and  temperature,  stress  jump  tests  must  be  conducted  (Follansbee,  Kocks  and 

Regazzoni  1985).  After  each  dynamic  test  with  stop  rings,  the  specimens  must  be 

reshaped  as  soon  as  possible  to  conduct  quasi-static  tests  without  any  room  temperature 

annealing.  Therefore  it  would  take  considerably  more  time  to  complete  tests  required  for 

MTS  model  compared  to  regular  tests  required  for  other  material  models.  Meanwhile  it 

has  more  than  fifteen  parameters,  which  makes  it  challenging  to  fit  all  parameters 

computationally  without  in  advance  assuming  values  in  some  of  the  parameters. 

Even  though  ZA  model  is  considerably  simpler  than  the  MTS  model,  it  is  not 

possible  to  model  the  aluminum-copper  alloys  due  to  its  hardening  term.  ZA  model  has  a 

1 

hardening  term  of  B0£2e~aT  ,  thus  flow  stress  is  proportional  to  square  root  of  the  strain; 
while  aluminum-cupper  alloys  show  saturation  behavior  after  some  point(Fig  3.1).  But 
hardening  parameter  can  be  modified  to  show  saturation  stress. 
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Figure  3.1.  Polycrystalline  aluminum  shows  hardening  behavior  akin  to  ZA  model’s 
hardening  parameter,  whereas  aluminum-copper  alloys  show  decreasing  hardening 
rate  with  increasing  strain  rate. 


If  original  Zerilli- Armstrong  model  is  used  to  model  the  Al-Cu-Mg-Mn-Ag  alloy,  it  is 
not  possible  to  capture  the  increased  rate  sensitivity  after  1000/s  strain  rate  in  the 
dynamic  regime  and  also  the  error  between  the  model  and  the  experimental  data  grows 
with  increasing  strain  which  indicates  a  different  strain  hardening  behavior  than  the  ZA 
model  (fig.  3.2). 

3.2  Former  Modifications  to  ZA  Model  by  Zerilli  and  Armstrong 

To  overcome  the  difficulty  of  modeling  materials  that  show  saturation  behavior, 
Zerilli  and  Armstrong  modified  the  ZA  model  to  include  dynamic  recovery  process 
(Zerrilli  and  Armstrong  1997).  The  mechanical  strength  of  obstacles  to  dislocation 
glide,  <7,  the  threshold  stress,  is  related  to  dislocation  density  : 

d  =  apbpz  (3-1) 


Therefore  dislocation  density  can  be  used  to  model  the  threshold  stress.  Heuristically,  the 
density  of  the  immobile  dislocations  may  be  related  to  the  strain  by: 
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Figure  3.2.  Original  Zerilli-Armstrong  Model  is  used  to  model  the  Al-Cu-Mg-Mn-Ag 
alloy.  Model  fails  to  capture  the  increased  strain-rate  sensitivity  in  the  dynamic 
regime. 


where  A  is  the  mean  free  path  for  immobilization  of  mobile  dislocations  and  co  is  the 
probability  for  remobilizing  or  annihilating  a  stopped  dislocation.  Eq.  3.2  can  be 
integrated  to  find  the  relation  between  strain  and  dislocation  density: 
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dp 

-j - =  ds 

JE-^P 

[  y~^ —  =  I  da  (3-3) 

Ib~“P 

1 

—  ln(  1  —  Abcop)  —  —e 


By  merging  Eq.3.3  with  Eq.  3.1,  an  expression  for  the  threshold  stress  can  be  found: 


—  (1  —  e-"£)  (3.4) 


With  this  new  definition  of  the  threshold  stress,  the  thermal  component  of  the  ZA  model 
becomes: 


where 


oth  =  Be  PT  + 
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_(1  -  e-"£) 

O) 
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—aT 


a  —  a0  —  ax  ln(f) 
P  =  Po  ~  Ei  ln(e) 


(3.5) 


(3.6) 


A  different  problem  with  the  ZA  model  in  modeling  aluminum-copper  alloys  is 
capturing  the  increased  strain  rate  sensitivity  of  these  alloys  in  the  dynamic  regime. 
Zerilli  and  Armstrong  again  modified  the  thermal  component  of  ZA  model  to  include  this 
strain  rate  sensitivity  effect  in  the  dynamic  region: 
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where  c0  is  a  constant  in  the  order  of  10  5  MPa/K.  In  the  quasi-static  regime, 


4  C0eT 

- term 

Oth 


will  be  very  small  and  negligible,  and  therefore  o*  will  be  almost  equal  to  oth.  However, 
as  4c°eT  term  gets  bigger  and  bigger  in  the  dynamic  region,  the  ZA  model  will  show  the 

°th 

desired  increased  rate  sensitivity  in  the  dynamic  region. 

3.3  Modifications  to  ZA  Model 

With  the  entire  modifications  ZA  model  become  (Modified  Zerilli-Armstrong 
Model): 
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where 


a  =  a0  —  a±  ln(f) 


(3.9) 


P=Po~Pi  ln(^) 

With  these  modifications  to  the  MZA  Model,  the  model  finally  can  capture  the 
increased  strain-rate  sensitivity  in  the  dynamic  regime.  As  Figure  3.3  illustrates,  the  MZA 
model  can  accurately  capture  the  experimental  data  in  quasi-static  and  dynamic  regimes 
at  room  temperature  and  220°  C.  Although  the  equations  seem  to  work  pretty  well,  the 
problem  is  with  the  elevated  temperatures.  Flow  stress  is  associated  to  the  temperature  in 
an  exponential  manner  which  gives  unrealistic  results  after  some  point  (Fig  3.4).  To 
accurately  capture  the  experimental  data  in  a  greater  temperature  range,  the  exponential 
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relation  between  the  flow  stress  and  the  temperature  must  be  altered.  Only  way  to  modify 
this  relation  is  to  make  one  of  the  constants  temperature  dependent.  Constants  a,  /. i ,  b  are 
more  or  less  constant  with  temperature,  and  the  term  B  cannot  be  assumed  to  decrease 
with  temperature  as  there  is  not  any  physical  basis. 
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Figure  3.3.  Al-Cu-Mg-Mn-Ag  alloy  modeled  with  Modified  Zerilli- Armstrong  Model 
which  shows  increased  strain-rate  sensitivity  in  the  dynamic  regime. 

Neither  of  the  two  physical  material  models  (MTS  and  ZA  models)  takes  account 
of  the  effect  of  the  diffusion  process  on  the  threshold  stress.  With  the  help  of  the 
diffusion,  dislocations  can  overcome  the  barriers  by  leaving  their  planes. 
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Figure  3.4.  Thermal  softening  of  the  flow  stress  is  shown.  After  500K,  model  behavior  is 
very  dissimilar  to  the  experimental  results. 

Mean  free  path  for  immobilization  of  mobile  dislocations  (A)  will  be  affected  by 
this  process  by  getting  larger  with  increasing  temperature.  Thus  thermal  component  can 
be  modified  as: 
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where  A(T)  is  a  material  specific  temperature  dependent  constant  which  decreases  with 
increasing  temperature.  Also  if  Be“^r  term  is  removed,  thermal  stress  becomes: 
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With  this  modification,  Be~^T  tenn  can  be  removed  since  A(T)e~aT  component  of  the 
new  thermal  stress  will  work  same  as  the  latter,  and  also  with  this  modification  threshold 
stress  becomes  analogous  to  threshold  stress  definition  in  MTS  model  (Eq.  1 .49  and  Eq. 
1.50),  where  threshold  stress  has  strain  dependent  and  independent  components. 
Furthermore  by  removing  this  term,  the  new  model  will  have  three  less  parameters 
(/?i,/?2,and  B)  without  any  loss  in  functionality  and  deviation  from  dislocation  physics. 
Considering  all  modifications,  Turkkan-Vural  Modified  Zerilli-Armstrong  (TVZA) 
model  turn  into: 
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where 


(3.13) 


After  these  modifications,  ZA  model  can  finally  capture  the  material  behavior  of  Al-Cu 
alloys  at  high  temperatures  and  A(T)  can  be  modified  to  capture  the  desired  thermal 
softening  behavior  of  other  different  materials(Fig.  3.5).  Similar  to  the  MZA  Model 
which  was  effective  at  capturing  the  experimental  data  at  relatively  low  temperatures, 
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TVZA  model  can  accurately  capture  the  experimental  data  at  room  temperature  and  220° 
C. 
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Figure  3.5.  TVZA  model  can  successfully  capture  the  softening  behavior  of  the  alloy 
even  at  temperatures  above  500  K. 
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Figure  3.6.  Al-Cu-Mg-Mn-Ag  alloy  modeled  with  Turkkan-Vural  Modified  Zerilli- 
Armstrong  Model  which  shows  increased  strain-rate  sensitivity  in  the  dynamic  regime. 
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CHAPTER  4 

MODELING  THE  EXPERIMENTAL  DATA 


4.1  Introduction 

Al-4.5%Cu,  Al-4.5%Cu-0.5%Mg,  and  Al-4.5%Cu-0.5%Mg-0.5%Mn  alloys  have 
only  room  temperature  data,  but  the  rest  of  the  materials  have  room  temperature  and 
elevated  temperature  experimental  data  in  quasi-static  and  dynamic  regimes.  Therefore 
the  following  alloys  which  have  complete  set  of  experimental  data  were  modeled  in  this 
chapter: 

•  [100]  and  [111]  Single  Crystals 

•  Polycrystalline  Al. 

•  Al.-0.1%Cu. 

•  Al-4 . 5  %Cu-0 . 5  %Mg-0 . 3  %Mn-0 .3  %  Ag 

As  discussed  in  Chapter  3,  original  Zerilli-Annstrong  Model  excels  in  capturing 
less  complex  alloys  such  as  Single  Crystal  Aluminum,  Polycrystalline  Aluminum  and 
Al.-0.1%Cu.  Therefore  these  materials  were  modeled  using  the  original  ZA  model  with 
only  increased  rate  sensitivity  modification  (Eq.  3.7). 
o  =  aa  +  o* 
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where 
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a  —  a0  —  ax  ln(d) 


(4.2) 


Al-4.5%Cu-0.5%Mg-0.3%Mn-0.3%Ag  alloy  was  modeled  using  the  Turkkan- 
Vural  Modified  Zerilli- Armstrong  Model.  Moreover,  the  listed  alloys  were  modeled 
using  Modified  Johnson-Cook  Model. 

4.2  [100]  and  [111]  Single  Crystals 

Single  Crystalline  Aluminum  has  only  room  temperature  data  and  no  further 
experiments  will  be  conducted  at  elevated  temperatures.  Therefore,  Single  Crystalline 
Aluminum  was  modeled  using  only  original  ZA  Model  (Fig.  4.1)  (Fig.  4.2). 


Single  Crystal  [111  Direction]  Single  Crystal  [111 

0.001/s  Strain  Rate  Direction]  3400/s  Strain  Rate 


Figure  4.1.  Single  Crystalline  [111]  direction  aluminum  modeled  with  original  ZA 
Model. 
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Figure  4.2.  Single  Crystalline  [111]  and  [011]  direction  aluminum  modeled  with  original 
ZA  Model. 

4.3  Polycrystalline  Aluminum 

Polycrystalline  Aluminum  have  complete  set  of  experiments  in  room  temperature, 
120°  C  ,and  220°  C.  Modified  Johnson-Cook  Model  and  Zerilli- Armstrong  Model  were 
used  to  model  the  experimental  data. 
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Figure  4.3.  Polycrystalline  aluminum  modeled  with  MJC  Model. 
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Figure  4.4.  Polycrystalline  aluminum  modeled  with  ZA  Model. 


4.4  Al-0.1%Cu  Alloy 

Al-0.1%Cu  alloy  have  complete  set  of  experiments  in  room  temperature,  120°  C 
,and  220°  C.  Modified  Johnson-Cook  Model  and  Zerilli- Armstrong  Model  were  used  to 
model  the  experimental  data  (Fig  4.5)(Fig.  4.6). 
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Figure  4.5.  Al-0.1%Cu  modeled  with  MJC  Model 
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Figure  4.6.  Al-0.1%Cu  modeled  with  ZA  Model. 


4.5  Al-4.5%Cu-0.5%Mg-0.3%Mn-0.3%Ag 

Al-4.5%Cu-0.5%Mg-0.3%Mn-0.3%Ag  alloy  have  complete  set  of  experiments  in 
room  temperature,  120°  C  ,and  220°  C.  Also  tests  were  done  at  elevated  temperatures  in 
0.001/s  strain  rate.  Modified  Johnson-Cook  Model  and  Turkkan-Vural  Modified  Zerilli- 
Armstrong  Model  were  used  to  model  the  experimental  data  (Fig  4.7)(Fig.  4.9).  Figure 
4.8  and  Figure  4.10  display  the  thennal  softening  behavior  captured  by  the  two  models. 
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Figure  4.7.  Al-4.5%Cu-0.5%Mg-0.3%Mn-0.3%Ag  modeled  with  MJC  Model. 
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Figure  4.8.  Al-4.5%Cu-0.5%Mg-0.3%Mn-0.3%Ag  thermal  softening  captured  with  MJC 
Model 
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Figure  4.9.  Al-4.5%Cu-0.5%Mg-0.3%Mn-0.3%Ag  modeled  with  TVZA  Model. 
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Figure  4.10.  Al-4.5%Cu-0.5%Mg-0.3%Mn-0.3%Ag  thermal  softening  captured  with 
TVZA  Model 
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CHAPTER  5 

ADIABATIC  SHEAR  BANDING 


5.1  Introduction 

As  discussed  at  section  1.4,  shear-compression  specimens  (SCS)  were  used  to 
study  adiabatic  shear  banding  on  different  alloys.  The  following  alloys  were  tested  with 
SCS: 

•  A12139-T8 

•  Al-4 . 5  %Cu-0 . 5  %Mg-0 . 3  %Mn-0 .3  %  Ag 

•  Al-4 . 5  %Cu-0 . 5  %Mg-0 . 3  %Mn 

5.2  Parameter  Fitting 

The  constants  k1  and  k2  (Eq.1.27  and  Eq.1.28)  are  found  by  comparing 
cylindrical  non-slotted  specimen  test  with  shear-compression  specimen  test,  both  of 
which  are  carried  out  at  similar  quasi-static  strain  rates. 


plastic  strain 


Figure  5.1.  Constants  /q  and  k2  are  found  by  comparing  quasi-static  SCS  and  cylindrical 
experimental  results. 


59 


Al-Cu-Mg-Mn-Ag  SCS  parameter  fitting  -  0.001  Is  Strain 


plastic  strain 


Figure  5.2.  Constants  k1  and  k2  are  found  by  comparing  quasi-static  SCS  and  cylindrical 
experimental  results. 


Al-Cu-Mg-Mn-Ag  SCS  parameter  fitting  -  0.001  Is  Strain 


Figure  5.3.  Constants  kx  and  k2  are  found  by  comparing  quasi-static  SCS  and  cylindrical 
experimental  results. 
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5.3  Dynamic  Testing 

Three  alloys  were  tested  using  SHPB  in  similar  strain-rates  (Fig  5.4)  and  Al- 
4.5%Cu-0.5%Mg-0.3%Mn  alloy  have  shear  localization  at  largest  strain  (0.45-0.5)  value 
among  other  alloys.  The  alloys  Al-4.5%Cu-0.5%Mg-0.3%Mn  0.3%Ag  and  A12139-T8 
show  shear  localization  around  0.35  to  0.4  strain.  But  these  tests  are  not  enough  to  assess 
shear  localization  behaviour  of  the  individiual  alloys,  and  therefore  stop-ring  tests  must 
be  done  to  investigate  microstructure  during  shear  localization  . 


0  0.1  0.2  0.3  0.4  0.5  0.6  0.7 


eq.  plastic  strain 


Figure  5.4.  SCS  results  of  the  three  tested  alloys.  Strain  rate  is  around  4000/s  for  all  tests. 

5.4  Stop-Ring  Testing 

To  study  how  shear  localization  develops  in  different  alloys,  dynamic  tests  need 
to  be  stopped  at  three  different  points;  first  at  the  beginning  of  the  deformation,  then  just 
before  the  shear  localization  initiates,  and  finally  after  the  localization  point  before  the 
fracture.  All  tests  except  the  ones  after  localization  point  are  complete  (Fig  5.5). 


eq.  true  stress[MPa] 
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Stop-Ring  Results 


Figure  5. 5. Stop-ring  SCS  results  of  the  three  tested  alloys.  Strain  rate  is  around  4000/s  for 
all  tests. 
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CHAPTER  6 
DISCUSSION 


6.1  Discussion 

Current  study  on  the  rate  and  temperature  dependent  mechanical  behavior  of  Al- 
Cu  material  system  found  out  that  the  Zerilli-Armstrong  (ZA)  model  is  not  sufficient  to 
model  even  the  simplest  aluminum  based  materials  without  increased  strain  rate 
modification  in  dynamic  regime.  After  this  modification  discussed  in  Section  3.2,  the 
constitutive  response  of  simple  materials  such  as  single  crystal  A1  (along  [001]  and  [111] 
directions),  polycrystalline  Al,  and  Al-0.1%Cu  can  be  modeled,  but  the  hardening  tenn  in 
ZA  model,  which  is  proportional  to  the  square  root  of  the  strain,  still  remains  to  be 
unrealistic  for  more  complex  alloys.  Even  with  the  Zerilli  and  Armstrong’s  alteration  of 
strain  hardening,  mechanical  behavior  of  more  complex  precipitation  hardened  Al-Cu 
alloys  cannot  be  captured  at  elevated  temperatures  due  to  exponential  relation  between 
flow  stress  and  temperature.  Current  study  proposes  further  modifications  to  modified  ZA 
Model  by  incorporating  the  effect  of  diffusion  and  precipitate  spacing  at  elevated 
temperatures  to  better  capture  the  thermal  softening  of  flow  stress  observed  in  current 
experiments. 

Another  objective  of  this  study  was  to  investigate  the  rate  and  temperature 
dependent  mechanical  behavior  of  Al-Cu  alloy  systems  and  to  this  end,  single  crystal  Al 
(along  [001]  and  [111]  directions),  polycrystalline  Al,  Al-0.1%Cu,  and  Al-Cu-Mg-Mn- 
Ag  materials  were  tested  and  characterized  at  room  and  elevated  temperatures.  MJC 
Model  and  the  newly  developed  TVZA  model  were  proposed  to  describe  rate  and 
temperature  dependent  constitutive  response  of  these  materials.  These  materials  show 
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similar  trends  in  their  material  response,  such  as  a  distinct  change  in  strain  rate  sensitivity 
in  quasi-static  and  dynamic  deformation  regimes.  Also  At  room  temperature  and  quasi¬ 
static  rates  (10‘3  s’1  to  10°  s’1),  flow  stress  is  almost  independent  of  the  stain  rate  while  in 
dynamic  regime  there  is  a  significant  increase  in  strain  rate  sensitivity.  However,  as  the 
temperature  is  elevated  above  room  temperature,  strain  rate  dependence  is  observed  even 
under  quasi-static  deformation  conditions.  On  the  other  hand,  in  dynamic  regime  (10"  s" 
to  104  s’1)  temperature  does  not  seem  to  influence  rate  sensitivity.  This  complex  behavior 
was  captured  with  phenomenological  MJC  Model  and  the  dislocation  physics  based 
TVZA  Model. 

Within  the  scope  of  this  project,  Al-Cu.,  Al.-Cu.-Mg.,  Al.-Cu.-Mg.-Mn.  alloys 
were  tested  at  room  temperature  and  due  to  lack  of  high  temperature  data,  these  alloys 
have  not  been  characterized  with  material  models. 

Last  aspect  of  this  study  was  to  assess  the  ballistic  properties  of  the  complex 
alloys  and  henceforth,  SCS  testing  was  proposed  to  study  the  tendency  of  the  alloys  to 
shear  localization.  SCS  results  revealed  that  Al.-Cu.-Mg.-Mn.,  Al.-Cu.-Mg. -Mn.-Ag.  , 
and  A12139-T8  alloys  show  shear  localization  after  similar  strain  values.  Furthennore 
stop  ring  tests  were  conducted  to  better  evaluate  shear  deformation  field  in  these  alloys. 
6.2  Recommendations  for  Future  Work 

The  following  is  recommended  for  future  work  for  the  characterization  of 
aluminum-copper  alloy  system: 

•  Complete  the  elevated  temperature  tests  for  the  rest  of  the  alloys  and  then, 
characterize  their  material  response  with  MJC  and  TVZA  Models. 
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•  Since  the  complex  materials  show  adiabatic  shear  banding  and  localization  at 
similar  strain  rates,  conduct  SCS  tests  at  elevated  temperatures  to  better  study 
the  ballistic  properties  of  the  complex  alloy  systems. 


APPENDIX  A 


MODELLING  PARAMETERS 
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Table  A.l.  Modified  Johnson-Cook  Model  parameters  for  different  materials 


Material 

Q 

o-o 

B0 

To 

T 

1  m 

P 

n  C± 

c2 

£01 

Et  k 

Poly.  A1 

2.9 

14.5 

140 

175 

780 

1 

0.45  0.0052 

0.1565 

0.001 

20  400 

Al-O.lCu 

3.9 

5.5 

112 

41 

670 

3 

0.36  0.0068 

0.2065 

0.001 

150  400 

Al-4.5Cu- 

0.5Mg-0.3Mn- 

0.3Ag 

4.5 

426 

234 

0 

640 

3.9 

0.2  0.0027 

0.07 

0.001 

200  400 

Table  A. 2.  Zerilli-Annstrong  Model  parameters  for  different  materials 

Material 

°a 

A 

Co 

a2 

Poly.  A1 

15 

424 

0.00000158489 

0.0044 

0.000225 

Al-O.lCu 

17 

390 

0.00000501187 

0.0039 

0.000175 
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Table  A. 3.  Turkkan-Vural  Modified  Zerilli-Armstrong  Model  parameters  for  different  materials 
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